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1. INTRODUCTION 

In the last few decades, our World has known a remarkable rise of the annual rate of population and 
urbanization which is estimated to 1.2%. As result, the rate of energy demand especially domestic electrical 
energy has grown up consequently over the 2009-2019 period to reach 23 430 TWh in 2019. In 2020, the 
world experienced a new tragic pandemic COVID-19 which turned our daily lives into disorder but despite 
all of that, this pandemic had a positive impact in the global primary energy consumption which fell by 4.5% 
that was considered as the largest decline since 1945. Also, the domestic electrical demand was declined for 
the first time since 2009 by 1.1% to reach 23 176 TWh as shown in Figure 1, according to enerdata statistics 
and statistical review of world energy 2021 [1], [2]. 

The global carbon emission from energy use was also estimated to have fallen by 6.3% since 2011 
as a positive result of the pandemic [2]. Staring from what’s happened last year, knowing the importance of 
reducing energy demand especially the non-renewable energy because of their negative environmental 
impact and the fact they are in limited supply. Also, since 2015 in the Paris Conference of Parties (COP21) 
there has been a huge interest and ambitions to decarbonize by most of countries which lead them to carry 
investigation in order to minimize the use of fossils fuels in electrical energy production and maximize the 
contribution of renewable energies such Wind, solar, geothermal. So, despite the disorder of 2020, renewable 
energy, continued to have a remarkably growth, wind and solar capacity increased by 238 GW last year 
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which is 50% higher than other last years. In the same way, the share of wind (173 TWh) and solar (148 
TWh) generation in the global power mix have known the biggest increase ever. Other way to increase and 
ameliorate the participation of Renewable energy in responding to the global electrical energy demand is 
encouraging each country to do its role. Indeed, in response to that our country Morocco decide to exceed 
objectives announced at COP21 by increasing the share of energies to more than 52% in the electricity mix 
by 2030 for that we have 50 renewable projects with a total installed capacity of around 4 GW that are 
already in operation. Morocco has more than 40 international companies that are active in a renewable power 
generation market. Since 2009, more than 12 billion dollars are invested in renewable energy. These projects 
have contributed to improve our energy independence by more than 8% since 2009, other projects are being 
developed or under construction and the multiple others will be financed over the next decade. In the same 
way Morocco encourage researchers to do research and studies in this fields thing that going to be subject of 
our article. 
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Figure 1. Electrical energy demand over 1990-2020 (TWh) [1] 


Wind is considered as one of the most renewable sources of energy. It is used to generate electricity 
from kinetic energy created by reason of the air motion; this generation is taking place using an electrical 
system called wind energy conversion system which transforms the wind energy into mechanical energy and 
then to an electrical one and that through many electrical and mechanical sub-systems that will be more 
developed in this article. The electrical energy generated by wind systems has known an enormous growth 
over the last years and its penetration into the power system have increased consequently so according to 
wind power statistics of the World Wind Energy Association, the installed wind capacity has reached 744 
GW by end of 2020 [2]. The Figure 2 shows the worldwide evolution of the total cumulative installed 
capacity through the 6 last years. 
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Figure 2. Total cumulative installed capacity over 2015-2020 [2] 
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To generate electricity from wind we have to pass by four essential sub-systems, described at Figure 3, 
so at first the wind hits the turbine’s blades and makes it rotate around a rotor that is connected to a gearbox 
which multiple the number of tours to an adaptive one and transmit it to the generator through the drive shaft, 
this generator converts the mechanical power into an electrical one that is injected to the grid through power 
converter. The combination of all these elements forms the wind energy conversion system (WECS) [3]. The 
total cumulative f power produced by the systems can vary depending on the size of the turbine and the 
length of its blades. The output depends on the rotor dimensions and the cube of the wind speed. So 
automatically, when wind speed increase, wind power generated increases. 


Wind : Multiplied 
>| Turbine Mechanical) Gearbox+ |Mechanical| Generator | Electrical | Power | Electrical Power 
Power Drive shaft Power Power Inverter to the Grid 
Figure 3. WECS blocs 


Since the power inverter plays a vital role in the WECS by connecting the power generator to the 
electrical grid in order to inject the electrical power generated and by converting DC power to AC and AC to 
DC to be suitable for the bidirectional power flow; There have been many developments and ameliorations in 
power electronic converters [4]. These converters are commonly used to regulate and adapted the electrical 
signal to the required form and parameters for a better integration in the grid and that’s by respecting the 
IEEE conditions for grid integration which consist on limiting of the total harmonic distortion (THD) in less 
the 5%. Indeed, power inverters become subject of several work and studies in order to attain the best 
topology which would make the penetration of the generated power easier and better. 

In 1980’s, as result of researchers’ investigations they have developed a new configuration of 
inverters called Multilevel inverter [5], [6]. A multi-level inverter is a power electronic system which 
provides an output alternating voltage level by using a number of DC voltage in input depending on the 
number of levels that will be obtained. The notion of multi-level inverter is a change made on the 
conventional two level one which gives an output wave form that varies between two values. So, the aim of 
multi-level inverter is to generate a smoother output waveform depends on the level, more level means that 
the output voltage will be smoother and near to a sinusoidal output. Multilevel inverters are very promising 
because of the many advantages they present as having a current with a reduced THD factor, lower stress of 
the electronic components due to decreased voltages, switching losses that are lower than those of 
conventional two-level inverter and also a smaller filter size is required, and lower electromagnetic 
interference (EMI), all of that make the multilevel inverter low-priced than the conventional one [7]. 

The purpose of this article is to present a comparative study and analyze of 3 different multilevel 
topologies of inverter that will be used in our wind energy conversion system and that’s in a way to prove 
that the current’s THD will be reduced by increasing the level of the multilevel inverter for a better 
integration in the grid also by improving and smoothing the output voltage waveform. Also, for a better 
penetration of the grid a control methodology of active and reactive power is proposed. The current paper is 
organized as follows. The introduction is presented in the first section. The section 2 presents a modeling of 
all WECS subsystems. Then, the analyze and model of two-level, three-level and five-level inverters are 
given in the section 3. After that we present in section 4 the control strategy used in our WECS. In the section 
5 simulations results in MATLAB/Simulink are presented and discussed. Last section concludes about the 
results obtained and gives an overall view of our further works. 


2. DYNAMIC MODELING OF WIND TURBINE SYSTEM USING A DFIG 

The aim of a wind energy conversion system is transforming wind kinetic energy into an electrical 
one. WECS is composed generally of a wind turbine transforming the wind energy into a mechanical one, a 
gearbox which adapts the rotation speed, a generator that transforms the mechanical power to an electrical 
one and power electronics in order to inject the generated power directly to the grid. In this section, we are 
presenting the model of each part of this system as in [8]. First, we are modeling the aerodynamic wind 
turbine that is composed of a horizontal three-blade and blades rotor which is linked to the gearbox using a 
low-speed shaft. Then we will present the model of the used generator which is a doubly fed asynchronous 
generator (DFIG) that we have selected for the many advantages it presents comparing to other generators. 
This one is directly connected to the grid using two power multilevel converters and for a better integration the 
use of controller is essential in our system; we are using the PI Controllers in both rotor side and grid side with 
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MPPT controller for a maximum power extraction and a pitch control. Our WECS is described in Figure 4 and 
power system data is summarized in the Table 1. Table 1 shows the parameters of our wind energy 
conversion system. 


Turbine 
—> 
p Gear-Box 690V20Hz Transformer 
wind 
—> 
_ 


Multi-level Multi-level 


Figure 4. Global scheme of WECS connected to the grid 


Table 1. Turbine, DFIG, and grid side parameters 


Parameters Values Parameters Values 
Rated power 1.5 MW Stator Resistance R, 2.65 mQ 
Rated wind speed 12 m/s Rotor Resistance R, 2.63 mQ 
Density of Air 1.225 kg/m? Stator Leakage Inductance Ly, 0.1687 mH 
Blade radius 31m Rotor Leakage Inductance L,, 0.1337 mH 
Gearbox ratio G 59 Magnetizing Inductance Lm 5.4749 mH 
Rated Power Pn 1.5 MW DC-Link Voltage Vac 1320 V 
Rated Voltage Un 690 V DC-Link Capacitor C 10028.7 uF 
Nominal Frequency 50 Hz Filter Resistance Rr 0.1838 Q 
Rated Rotor Speed 1750 tr/min Filter Inductance Lp 0.61187 mH 
Number of Pole Pairs 2 


2.1. Wind turbine dynamic modeling 

The aerodynamic model of the wind turbine system has been studied and presented in several works 
[9], [10]. This part aims to transform the incoming air into a mechanical power using the aerodynamic 
equation below and to transmit it to the generator through the gearbox, the Figure 5 present the 3D 
characteristic of aerodynamic power: 


Paero = Cy, B) Py = Cp, B) - P Awt w? (1) 


P, is the wind’s kinetic power; Cp is the power coefficient depending on the blades pitch angle 6 and the 
turbine tip speed ratio A; Awt, p and w, are the turbine blades area, the air density and the wind speed. 
The tip speed ratio is given by (2): 


R Qeur į 


janmi (2) 


We i 


Neur ; is the speed of turbine, and R is the radius of turbine blades. 
Mechanical angular speed @me, and Gearbox torque T, are expressed as follow, where G is the 
gearbox coefficient. 


r = a) 


Wmec = G Wtur (4) 


The mechanical equation of our system is: 


dwmec 


J dt = Ty = Tem z f Wmec (5) 
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Tem and f are the DFIG electromagnetic torque and the friction coefficient. J represents the total inertia that 
can be expressed as (6): 


J= ur + Joric (6) 


G2 


Jpricg>Jtur are the inertias of the generator and the turbine. 
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Figure 5. 3D characteristic of aerodynamic power 


2.2. DFIG dynamic model 

To convert the mechanical power Pm produced by the turbine to an electrical power, a generator 
system is required. Several types of generators exists but for this work we are using the doubly fed induction 
generator DFIG because it is the most used system for variable speed WECS as studied in [11]—[14]. In this 
part, we are presenting DFIG equations in park reference dq linked to a rotating field. For simplifying those 
equations, we use a Park reference linked to the rotating field and the hypothesis of a stator flux is constant 
and oriented along the axis d and stator resistance as negligible. The stator and rotor generator dq equivalent 
model is presented in Figure 6. 


Figure 6. DFIG equivalent model in dq 


The DFIG model is represented by the following equations: 


Vsa R; 0 0 0 isa Psa 0 —Ws 0 0 Psa 
Vsq a3 0 R, 0 0 icq d Poq + Ws 0 0 0 Poq (7) 
Vra 0 0 R, Offing]* al%al7}0 0 0 -sa| Pra 
Urq 0 0 0 Rly Yal LO 0 so, 0 I%, 


i, v and are the current, the voltage and the DFIG field. L and R are the inductance and the resistance; s 
and r indicate the side stator and the side rotor of generator and w, is the electrical angular speed of the 
DFIG. Stator and rotor magnetic fields equation are: 


Reduction of harmonics emission of a WECS in the electrical ... (Maha Annoukoubi) 


1524 O ISSN: 2088-8694 


Wal T-L, 0O Lm O0Jfia 
Peq = 0 =L; 0 Lm isq (8) 
Pra 7 Lm 0 =L, 0 ira 
Woa O Im 0 -Ld liq 


Ls; L, and L,, are the stator inductance; the rotor inductance and the mutual inductance. Equations of 
electrical active and reactive powers delivered by the DFIG are expressed as: 


oleh alle | a 


P, 23 Pa a 
Al i dl Vra —Vra la (10) 


2.3. Mode of the grid side RL filter 
WECS inject the generated power to the grid through an RL filter, the model of this filter is 
expressed as follow: 


K 0jJafi —R Wel k 
hdo -idal Eat edil bu an 
94 hh tt ggd lL Weblp hiy ul 
where Vra and vf, are the output voltage of the filter in d-q ; i¢g and i¢g are the filter current; the voltages 
Vig and Vig are inverter’s output voltage in d-q frame; Ry and Ly are the resistance and the inductance of the 


filter. 
The resulted active and reactive powers (Pwt ,Qwz) that will be injected to the grid are expressed by: 


on] = lal+ H (12) 


3. MODEL OF MULTILEVEL INVERTERS AND THEIR IMPACT IN TERMS OF 
HARMONICS 

Recently, multilevel inverters have known an important development especially in terms of grid’s 
connection with WECS. A multilevel inverter is a power electronic system that enables obtaining a specific 
alternative output voltage from a DC input voltage; it is a new concept that brings new modification and 
amelioration to the conventional two-level inverter. The aim of multi-level inverter is to generate a smoother 
output waveform depends on the level, more level means that the output voltage will be smoother and near to 
a sinusoidal output. Multilevel inverters are very promising because of the many advantages presented at 
[15], [16] in terms of improving power quality and higher voltage capability; providing a current with a 
reduced THD factor; low THD; lower stress of the electronic components due to decreased voltages; 
switching losses that are lower than those of conventional two-level inverter; smaller filter size is required; 
lower EMI; and low-priced than the conventional one. 


3.1. Multilevel topologies 

Various topologies for multilevel inverters have been proposed since the development and 
improvement of the power semiconductors as the insulated gate bipolar transistor (IGBT) and metal oxide 
semiconductor field effect transistor (MOSFET); Common ones are diode clamped, flying capacitor or multi- 
cell, cascaded H-bridge, and modified H-bridge multilevel and they are all differentiated by their topologies 

and control methods as presented in Figure 7 [17]: 

— The neutral point clamped (NPC) inverter which was the first multilevel topology that was proposed by 
Nabae [5]. They are using clamping diodes for guaranteeing the proper voltage sharing across the power 
switches and as result, it decreased the THD factor in output voltage and current and smooth their 
waveforms. Despite the advantage of this topology but it is limited only in three levels and that due to the 
arrangement complexity of power semiconductors in higher levels. 

— Flying capacitor is a developed topology of NPC that replaced the use of diode by capacitor [18]. This 
structure gives more flexibility in waveform synthesis and balancing voltage. Redundancy in the 
switching states is available by using flying capacitors instead of clamping diodes. And it could be 
utilized to control the regulation of the capacitor voltages. 
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— H-bridge structures [19], [20], which has known a high interest due to the greater demand of medium- 
voltage high power inverters. The H-bridge cascaded inverter uses series of single-phase full-bridge 
inverters to construct multilevel phase legs wit separate dc sources. 


Multilevel Inverters 


Common DC Sources Separate DC 
Diode Clamped NPC Flying Capacitor H-bridge 


Figure 7. Multilevel inverters topologies 


3.2. Harmonic impact 

Harmonic emissions are identified as a main power quality problem for WECS, that’s why it is 
essential to understand harmonic behavior of WECS in order to control their effect on the electrical grids to 
assure a good connection. In the power system, the definition of a harmonic can be stated as follows: a 
periodic sinusoidal component having a frequency integer multiple of the frequency of the fundamental [21]. 
So, for a system, the fundamental frequency is fo and the harmonic h”! frequency is /fo. Harmonics are often 
used to define the distortion of the sinusoidal signal associated with current or voltage of different amplitudes 
and frequencies which results a degradation of the power factor. 

The concept of harmonic was introduced at the beginning of the 19th century by Joseph Fourier who 
demonstrated that any periodic non-sinusoidal signal can be represented by a sum or series of sinusoids of 
discrete frequencies, multiple. 


f= > + Ez (ancos(not) + b,sin(nwt)) (13) 


where Qo, an and b, are defined as: 


dg => J, f (wt)d(wt) (14) 
an = =." Ff (wt)cos(nwt)d(wt) (15) 
b, = Z J, f (wt)sin(nwt)d(wt) (16) 


The first component of rank 1 (h = 1) is called: fundamental component. For systems connected to a 
stabilized electrical grid, the frequency of the fundamental component is considered as fixed, 50 Hz in 
Europe and 60 Hz in the United States. The rest of the components of the series FOURIER are called: 
harmonics of order h, where h designates the number of the component, order 2 corresponds to the second 
term of the FOURIER series which will have a frequency double of the fundamental. Figure 8 gives an 
example of a current containing a third harmonic, i.e., a current which contains a rank 3 component high in 
amplitude. In three-phase electrical grid, the main harmonic components are of ranks 5, 7, 11 and 13. 


vy 
fundamental signal 
Distorted signal 


Harmonics 


T(s) 


Figure 8. Representation of the effect of a harmonic 
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The presence of current or voltage harmonics leads to negative effects on the distribution grid as 
studied in [22], such as: 

— Heating of conductors, cables, capacitors and machines due to additional copper and iron losses. 

— Electronic equipment may be sensitive to the voltage distortion supplying it, due to higher voltage peaks, 
unexpected zero-crossing, and affectation to protection circuits. 

— Reduced service life of components and equipment under continuous distorted supply voltage. 

— Resonance phenomena. The resonant frequencies of the circuits formed by the transformer inductors and 
the capacitances of the cables are normally quite high, but these can coincide with the frequency of a 
harmonic. In this case, there will be a significant amplification which can destroy the equipment 
connected to the grid. 

— Interference induced on communication lines, particularly electromagnetic radiation. 


3.3. Total harmonic distortion 

It is defined as a ratio between the r.m.s. value of all the harmonics and the r.m.s. of the fundamental 
frequency. It is usually given as a percentage [%] (by multiplying previous result per 100) [23]. Current 
THD: according to the definition above, the THD for current will be, 


THD, = {EN_,(#) a7) 


h 


Voltage THD (THDv): similarly, the THD of voltage is expressed as, 


THD, = [DX ey (18) 


n=2 Uy 


In order to respect the IEEE conditions for grid integration which consist on limiting of the THD in less than 
5%. The use of multilevel inverters seems to be a good solution to produce a stepped output phase voltage 
with a refined harmonic profile compared to a two-level inverter. 


3.4. Model of the two-level inverter 

A two-level inverter is a DC/AC converter that produces two different output voltages + Vdc/2 and 
— Vdc/2 from an input voltage Vdc and in order to create an AC voltage, these two generated voltages are 
continually switched. For switching we used a pulse width modulation (PWM). This method of generating 
AC output from a DC one can be effective but it has few drawbacks in creating harmonic distortions in the 
output voltage, it also has a high dv/dt comparing to the one providing by the multilevel inverter [24]. For 
this work we are using a three-phase two-level inverter that is composed of 3 arms each one consists of two 
switches controlled using a PWM Signal presented in the Figure 9. 


Figure 9. Configuration of a two-level inverter 


The two-level converter function matrix is presented as: 


1 —1 -1) [Sa 
-1 1 -1))Sp (19) 
-1 -1 1141S, 

In MATLAB/Simulink we modeled and simulated the two-level inverter as a part from the WECS 
in order to observe the output waveform and analyze the THD before using it with WECS. Figure 10(a) 
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present a single-phase two-level output voltage obtained from a DC voltage that is equal to 1320 V. 
Figure 10(b) is a zoom in the output voltage. THD result of the FFT analyze will be presented in the section 
of Simulation and results. 


o 
3 


800 


600 


Two-level voltage converter (V) 
o 
= — j 
Two-level voltage converter (V) 


1 L L L L L L L L - - 
0.1 0.2 0.3 0.4 0.5 0.6 07 0.8 0.9 1 0.52 0.54 0.56 0.58 0.6 0.62 0.64 0.66 O68 0.7 
Time (s) Time (s) 


(a) (b) 


Figure 10. Simulation results of two-level inverter, (a) Output voltage and (b) Zoom in the output voltage 


3.5. Model of the three-level inverter 

The configuration of the NPC three-level inverter is presented in the Figure 11, it is composed two 
capacitor that divide the Vdc input voltage in order to converter it into a three level AC output voltage and 
three arms each one is composed of four IGBT switches, two diodes as presented in [25]. Switching 
sequences for the three-level inverter are generated through PWM controlled by a control system presented in 
the next section. 


Figure 11. Configuration of three levels inverter 


The three-level inverter function matrix is: 


Vo4 P 1 —1 -1 Sar 2 1 -= —1 Saz 
Va| = 5Ucr|-1 1 —1||Sa|+zUa2|-1 1-11] So2 (20) 
Vs3 =1 =1 1 Ser =L =l 1 Sez 


Si are the Switching sequences applied by the PWM. In MATLAB/ Simulink we modeled and 
simulated the Three level inverter in order to observe the output waveform and analyze the THD before using 
it with WECS. Figure 12(a) present a single phase three level output voltage obtained from a DC voltage. 
Figure 12(b) presents a zoom in the results. THD result of the FFT analyze will be presented in the section of 
simulation and results. 
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Figure 12. Simulation results of three-level inverter, (a) Output voltage and (b) Zoom in the output voltage 


3.6. Model of the five-level multilevel inverter 

The five-level inverter used is an H-Bridge structure cascaded composed of three arms in cascaded 
each one is composed of two capacitors that divide the input Vdc in two and two H-bridge inverter each one 
of them has 4 IGBT. The model of a three-phases five level cascade H-Bridge inverter is presented in 
Figure 13 and also developed in [19]—[26] and the switching sequences of switches is presented in Table 2 . 


sik 4 Sul 4. sail 4. Su Ka sk A Soil 4. 


Ak KI | sk sa | sk sK 
: Vaj * Vaj * Vss 
sk 4 sulk 4 sik 4 silk A sik 4 sulk 4 
c= c G 
sal A sul as S: A a sal a Sal a sul > 


Figure 13. Model of three-phase five-level cascade H-bridge 


Table 2. Switching state 


Sı S2 S2 S, Sy! S2 53! Sy! Vout 
1 1 0 0 0 0 1 1 0 

1 0 0 1 1 0 0 1 0 

1 0 0 1 0 1 1 0 0 
0 1 1 0 1 0 0 1 0 
0 1 1 0 0 1 1 0 0 
0 0 1 1 1 1 0 0 0 

1 1 0 0 1 0 0 1 Vac 
1 1 0 0 0 1 1 0 Vac 
1 0 0 1 1 1 0 0 Vac 
0 1 1 0 1 1 0 0 Vac 
1 1 0 0 1 1 0 0 Woe 
0 0 1 1 1 0 0 1 Vac 
0 0 1 1 0 1 1 0 Vac 
1 0 0 1 0 0 1 1 Vac 
0 1 1 0 0 0 1 1 Vac 
0 0 1 1 0 0 1 1 -2V4 


In MATLAB/Simulink we modeled and simulated the five-level inverter in order to observe the 
output waveform and analyze the THD before using it with WECS. Figure 14(a) presents a single phase five 
level output voltage obtained from a DC voltage. Figure 14(b) gives a zoom for this result. THD result of the 
FFT analyze will be presented in the section of simulation and results. 
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Figure 14. Simulation results of five-level inverter, (a) Output voltage and (b) Zoom in the output voltage 


4. 


WECS CONTROL STRATEGIE OF ACTIVE AND REACTIVE POWER 

In this section, we present the control strategies used in our WECS in order to improve the quality of 
produced electrical energy for a better integration into the electrical grid. The global scheme of WECS and 
controllers is presented in the Figure 15. 
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Figure 15. Control scheme of WECS 


4.1. Maximum power point tracking algorithm 

In order to extract maximum power possible from the kinetic energy of wind, a maximum power 
point tracking algorithm has been developed. For our WECS we used MPPT controller in order to track a 
specific speed for extracting maximum power from the wind despite the wind speed variation. The principle 
of this controller is to keep the tip speed ratio equal to an optimum value in order to maximize the power 
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4.2. Pitch angle control 

Other way for extracting maximum power is managing the orientation of blade by regulating the 
pitch angle ß of blades. The blade orientation system consists on positioning the blades at a reference B’°S 
angle which depends on the direction of wind, this rotation is done by a servo motor. This system enables 
turbine blades to rotate using the Pitch control [27], [28]. The pitch angle 87°’ depends upon the mechanical 
power control P,. Pitch Angle f equation can be expressed as, 


B= pret (23) 


TF S 
s is a Laplace Operator; Tg is the constant time of the pitch servo. 


4.3. Control strategy of the rotor-side inverter 

For a better integration of the generated power into the grid a control of the used inverters is 
necessary, so in order to control the rotor-side inverter, we used the Indirect Stator Field Oriented Control 
(ISFOC), the principle of this control is orienting the stator field vector Y, along the axis d [29]. So 
that: Psg = 0; Ysa = Ys. As known, in high and medium power machines the stator resistance Rs can be 
considered as negligible, and considering that the grid is stable, so that the stator field ¥, is constant. 

DFIG model can be expressed as, 


|| El 
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The rotor angle expression is, 
0, = J (wg F Omec) dt (26) 


As we are using the ISFOC, the compensation and the coupling terms could be expressed as, 
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Electromagnetic torque equation is, 
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From the previous equations, we obtain the reference currents te f and ie for regulating the produced active 
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4.4. Control strategy of the grid-side inverter 

In order to control the grid side inverter, we used the voltage-oriented control (VOC) because and 
according to [30], [31] it is considered one of the efficient controller for the grid side thanks to its lower 


energy losses and current distortion. The VOC used a Vdc control loop and an inner current control loop in 
order to have a fast dynamic response, so, Vgqg = 0 and vg, = Vg. The equations of P; and Qp are, 


Pr = = Ugipg (33) 


3 ; 
Qr = 5 Vala (34) 


The filter currents are: 


d iga = =k, lgd 0 0 0 
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Using VOC we obtain a compensation and the coupling that are expressed as: 


CEE alle = 


lgq 


š ref - 
Direct reference current ae is: 


fa =n (37) 


For this WECS, we are using a PI controller for a separate control of the injected electromechanical 
torque and flux of the DFIG. PI controller is also used in order to ensure a follow of reactive and active 
powers of the DFIG to their references. In Figure 16 we present the block diagram of the closed loop for the 
PI control strategy. The general PI controller expression is: 
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Figure 16. PI Controller loop 


5. SIMULATION RESULTS AND DISCUSSION 

For evaluating and testing the performance of our WECS, equipped with a pitch controller, an 
MPPT control strategy, and two multilevel inverters connected to the rotor and to the grid and controlled by 
an ISFOC for the rotor side and a VOC in the grid side as shown in Figure 15. We carried out multiple 
simulations in MATLAB/Simulink environment to verify the performances of each multilevel separately, 
using a variable wind speed profile, which varies between 8m/s and 12m/s as shown in Figure 17. The 
studied WECS, DFIG and grid parameters are listed in Tables 1, 2 and 3. Depending on the result of the 
WECS simulation presented in Figures 18(a) and Figure 18(b), we can conclude that our system has two 
functioning mode the first one is the MPPT used for wind speed that is less than 12 m/s. So, when the wind 
speed exceeds the nominal speed, the aerodynamique power is limited using the pitch control regulation and 
so that the power coefficient and tip speed ratio decreases. 
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Figure 18. Simulation results of power coefficient and it’s parameters, (a) Power coefficient and (b) Speed 
ratio and pitch angle 


Figure 19(a) show that active power follows its reference as explained theoretically in the previous 
section. Figure 19(b) shows the same things for reactive power. Figure 20 present the output voltage Vdc of 
the rotor side inverter this result shows that it’s equal to its reference Vac rer = 1000 V. We also presented the 
injected grid voltage in blue and stator current in orange in the Figure 21, the current and voltage are not in 
phase and as we conclude, at the start-up time the current is equal to zero and then start to increase and varies 
depending on the wind profile 
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Figure 19. Simulation results of active and reactive power of the DFIG: (a) active power and (b) reactive 
power 
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This first part of simulation was only dedicated for the control of active and reactive powers of our 
WECS based on DFIG, a 5 level cascaded H-bridge multilevel inverter and a VOC and ISFOC controllers 
using a PI controller. The main objective of these controls is to extract the maximum power possible from the 
available wind speed in order to inject it safely into the grid. The obtained resulted of simulation by 
MATLAB/SIMULINK proved the efficacy of those controls strategies of regulating the active and reactive 
powers of our WECS. In this part we are going to highlights the performance of the 3 proposed multilevel 
inverter topologies studied in section 3 used for our WECS and that’s by elaborating a comparative study of 
the THD factor resulted from an FFT analyze done for the current that will be injected to the grid i, , the 
stator current is and the filter output current if . In Figure 22(a) we present the THD result using a two level 
inverter of the currents if, is in Figure 22(b) and i, in Figure 22(c). 
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Figure 22. FFT analyze result of currents of WECS using two-level inverter, (a) ir current, (b) is current, and 
(C) ig current 
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Figures 23(a), 23(b), and 23(c) presents the THD of if, is and i, using the three level inverter studied 
in section 3, and as we can see the THD factor has been reduced by 1.62% for current if, 0.67% for current is 
and 0.41% for current ig in comparison with the two level one. Figure 24(a), 24(b) and 24(c) presents the 
THD of if, is and ig using the cascaded H-bridge 5 level inverter studied in section 3, and as we can see the 
THD factor has been reduced by 0.02% for current if, 0.49% for current is and 2.59% for current i, in 
comparison with the three level one and by 1.65% or current if, 1.16% for current is and 3% for current i, in 
comparison with the result of the two-level inverter. We can conclude that the use of multilevel inverter in 
addition to the amelioration of the output voltage form studied in section 3 it decreased significantly the THD 
factor of the current that will be injected to the grid which is less than 5% and makes our WECS satisfy the 
IEEE 519 standards. Table 3 resumes THD results. 


Table3. THD results conclusion 
THD 1,(%) 1(%) _1,(%) 
2level 5.37 1.65 3.85 
3 level 3.75 0.98 3.44 
Slevel 3.72 0.49 0.85 
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Figure 23. FFT analyze result of currents of WECS using three-level inverter, (a) iy, (b) is, and (c) ig 
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Figure 24. FFT analyze result of currents of WECS using 5 level inverter: (a) is, (b) is, and (c) ig. 
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6. CONCLUSION 

This paper presents a comparative analyze and study of the performances of using multilevel 
inverter for WECS application. The theoretical study has been confirmed by an FFT analyze simulation done 
in MATLAB/Simulink which proved that THD factor of the current that will be injected to the grid have 
been reduced significantly using a 5 level inverter instead of three and two level ones. Moreover, the use of 
multilevel inverter for WECS instead of the conventional one present more advantages as a smoother output 
waveform which enable us to minimize the size of the filter or even eliminating it for higher levels also the 
dv/dt stress and EMI are reduced so that the dimensions of semiconductors used and consequently the cost of 
the inverter. By this theoretical and simulation results we aim to compare it with a practical one in order to 
confirm our study. 
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